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ABSTRACT 

We  state  upper  and  lower  bound  formulas  for  the  torsional  stiffness  of 
shafts  of  varying  circular  cross  section,  in  accordance  with  the  classical 
Michell  formulation  of  this  problem,  through  use  of  the  principles  of  minimum 
potential  and  complementary  energy.  The  general  results  are  used  to  obtain 
explicit  first-approximation  bounds  which,  for  the  limiting  case  of  the 
cylindrical  shaft,  reproduce  the  known  elementary  exact  results.  It  is  con- 
jectured that  the  first-approximation  lower  bound  is  significantly  closer  to 
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ON  BOUNDS  FOR  THE  TIORSIONAL  STIFFNESS  OF  SHAFTS  OF 
VARYING  CIRCULAR  CROSS  SECTION 


E.  Reissner 


We  consider,  in  terms  of  cylindrical  coordinates  r and  z,  a 

homogeneous  isotropic  linear  elastic  solid  bounded  by  plane  surfaces  z = 0 

and  z = L,  and  by  a surface  of  revolution  r = r (z)  for  0 « z £ L.  We  assume 
' o 

that  the  surface  portion  z = 0 is  held  fixed,  and  that  the  boundary  conditions 
for  the  surface  portion  z = L are  mixed  conditions,  of  prescribed  circumfer- 
ential displacement  v = ©r,  an^  of  vanishing  normal  stress  a and  shear 

stress  T . We  further  assume  that  the  surface  portion  r = r (z)  is  traction 
rz  o 

free. 

The  solution  of  this  problem  can  be  obtained,  as  first  shown  by  J.  H. 

Michell,  by  assuming  the  vanishing  of  all  radial  and  axial  displacements  and 

of  all  stresses  except  T a and  T With  this  there  remain  the  stress  strain 

r ro  z o 

relations  T Q = GT  and  T a = Gy  with  strain  displacement  relations 


yz6  =\ 


yr6  = r<r’lv>,r 


where  v is  independent  of  9,  and  with  stress  stress-function  relations 


T = r"3* 
zo  , r 


T,e*-'  *. 


where  & is  independent  of  9. 

While  the  boundary  conditions  for  z = 0 and  z = L are  now  directly 

given  in  terms  of  v,  the  corresponding  conditions  of  no  tractions  for 

r = r (z)  are  readily  shown  to  be  equivalent  to  the  one  condition  = const, 
o 

In  what  follows  we  are  interested,  in  particular,  in  the  values  of  the 

stiffness  coefficient  K in  the  torque -twist  relation  T = K©  where 

,r 

T = 2rrJ  °tz9  r* dr*  (3) 


1 


Specifically,  our  object  is  to  obtain  upper  and  lower  bounds  for  K,  through 
use  of  both  the  principles  of  minimum  potential  and  minimum  complementary 
energy,  analogous  to  what  has  recently  been  done  for  the  problem  of  the  end 
loaded  cantilever  beam,  treated  within  the  framework  of  the  theory  of  plane 
stress  [l]. 


The  appropriate  basic  inequalities  for  the  special  case  under  consid- 
eration are  readily  shown  to  be  of  the  form 
,.L  r 


(lL  r 

T©  - 2ir\  ° B rdrdz  s |T0  ^ 2tt\  1 ° A rdrdz. 

Vo 


In  this  we  have 
A 


= lG(?ze  S=2b(?z9  +frd)  • 


(4) 


(5) 


where  y and  V o are  given  in  terms  of  a function  v in  accordance  with 
z0  ru 

(11.  and  T and  T are  given  in  terms  of  a function  'I'  in  accordance 
' ' z 0 rb 

with  (2).  The  function  v must  be  differentiable  and  satisfy  the  boundary 
conditions  v(r,0)  = 0,  v(r,  L)  =©r.  The  function  'I'  must  also  be  differ- 
entiable  and  satisfy  the  boundary  condition  ^[r^zj.z]  = const,  for  0 ^ z 5 L. 
Furthermore,  v and  'I'  must  be  such  that  the  double  integrals  in  (4)  do  in 
fact  exist.  We  note  that  the  equal  signs  in  (4)  apply  in  the  event  that  v and 
>£,  respectively,  are  the  actual  solution  functions  v and  'I'  of  the  Michell 
boundary  value  problem. 


First  - Approximation  Upper  Bound.  We  assume 

v = r V(z)  (6) 

giving  = 0 and  y = rv'(z)  and  therewith 

L r _ 

II  = 2ir\  | A rdrdz  = \vG\  [V^z)]3  r4  dz.  (7) 

d Vo  *0  ° 

Wc  determine  V(z)  by  setting  = 0,  subject  to  the  constraint 

conditions  V(0)  = 0,  V(L)  =Q,and  obtain 


2 


V (z)  =0 


I r"‘ 

Jo  ° 

jV 

J „ O 


and  therewith 


4 I r 4 dz 
J0  ° 


and  then,  from  (4), 


U1  L _ o 

2 | r 4 dz 

Jo  ° 

We  note  that  is  the  approximate  value  of  K which  follows  by  an 

elementary  combination  of  the  solutions  for  a large  number  of  short  shafts 
of  appropriate  differing  uniform  circular  cross  section. 

First-Approximation  Lower  Bound.  We  assume 

* = c(r / r )4  . (11) 

o 

which  satisfies  ^(r  , z)  = c,  is  exact  for  the  uniform  diameter  case,  and  gives 

o 


-a  , 2 


4c  r r,  T a = 4c r r ' r 


Therewith 


n = T©  - 2ff  ! c 
3 Jrt«L 

0 0 


B rdrdz  = 2ttc©  - 


Ittc  j*  r 2 / a dz 

-77-  1 1 + 7 r J -j- . ( . 

G J „ 3 o r 

0 o 


We  determine  the  smallest  value  of  II  by  setting  ^11  />c  = 0.  This  gives 


4j  [i+ 1 v;f i C * 


• ,14» 

4l0ti  + 3 


3 


and  then,  from  (4), 
K 


rrG 

L1  2PL  [1  + (2/3)(r')2]  r‘4  dz 
0 ° 


(15) 


We  note  that  K = KTT  = K , as  it  must  be,  for  the  case  of  the  uniform 
LI  U 1 o 

circular  cylinder  for  which  r^  = 0,  throughout.  For  other  cases  we  expect, 

on  the  basis  of  experience  with  other  problems,  that  K will  be  a better 

i i 

approximation  than  K to  the  exact  value  of  K. 

Results  for  conical  shafts.  Setting  r^  = a + z tana,  where  a and 
a are  constants  we  find,  from  (10)  and  (15), 

K. 


K 


LI 

U1 


1 


1 + (2/  3)  tan2  a 


(16) 


showing  a significant  difference  between  K and  K for  sufficiently  large 

1~j  1 U i 

values  of  a.  For  example,  when  a - 30°  then  K /K  = 9/  11  = .82,  and 

*—j  1 U 1 

when  a = 45°  then  = .6.  We  note  that  the  known  closed-form 

solution  for  conical  shafts,  [2],  is  not  an  exact  solution  of  the  problem  as  stated 
here  inasmuch  as  this  closed-form  solution  implies  rigid  rotations  of 
spherical  end  surfaces,  rather  than  of  plane  end  surfaces. 

Higher-Approximation  Bounds.  It  appears  to  be  easier  to  derive 
improved  lower  bounds  of  relatively  attractive  appearance  than  to  do  the 
same  for  the  problem  of  upper  bounds. 

In  order  to  obtain  a sequence  of  lower  bounds  K for  N = 1,2... 


we  may  set 


$ = (r/r  )4[c  + c (r/r  )2  + ...  + c (r/r  )2N_21 

o 1 2 o No 


( IV) 


and  then  determine  values  of  the  coefficients  c from  the  N simultaneous 


linear  equations  an  lb c =0,  n = 1,2 

s " 

In  order  to  o 
generalization  of  (6) 


• N. 


In  order  to  obtain  improved  upper  bounds  we  may  set,  in 


4 


-w  2N  - 1 

v = 0 [ rVj.  (z)  + r3  V3  (z)  + + r V^z)!  * (18) 

In  this  the  functions  V have  to  satisfy  the  constraint  boundary  conditions 

^^(0)  =0,  Vj  (L)  = 1,  Vg(L)=0,...,V  (L)  = 0,  with  the  variational 

equation  611^=  0 then  leading  to  a system  of  simultaneous  linear  second- 

order  differential  equations  for  the  functions  V . 

n 

To  illustrate  the  nature  of  the  problem  we  consider  the  case 
v = 0(rV1  + r3V0)  where  we  will  assume,  additionally,  that  V,  concides  with 
the  solution  function  V in  our  first-approximation  calculation.  We  now 
obtain,  in  place  of  equation  (7) 

^d  = 4Ge3J'^[rS(Vi)2  +1  roV«  + 2 ro(Vz,)2  +T  roV^V^  dz  ' (19) 

and  with  this,  as  differential  equation  for  V3, 


<roV=>'  -T  r6oV=  =-f(,oV'>1'  - 


3f  I— ‘dz 

Jo  ° 


The  associated  minimum  value  of  Tld  follows  from  (19),  with  the  help  of  (20) 
and  (8),  in  the  form 


~ Ges  r.  4 „L  , t 

n , = ; 1 - ~ f r r V2^z]  • 

d 4 fL  r -4  dz  1 3 0 ° ° J 


and  then  from  (4),  and  with  K as  in  (10), 

o 


KU2  = Ko  [*  ‘I  Vo^H  * <22 

Equation  (22)  may  be  written,  with  the  help  of  (20)  and  upon  making  use  of 
the  boundary  conditions  for  Vs,  in  the  alternate  form 

KU2  = K„  {1-!f0L('e„Vl+TTr;iV^]dz}  . (22- 


ti 


which  makes  it  evident  that  Ky  is  in  fact  a better  upper  bound  for  K than  the 
first  approximation  bound  K^. 

As  regards  the  determination  of  the  function  V2  we  note  the  possibility 

of  an  explicit  solution,  as  a combination  of  powers  of  r , for  the  case  r = a 

o o 

+ z tana.  The  resulting  ratio  KU2  /Ku  j is  a less  simple  expression  than  the 
ratio  K^j/K  *n  061  and  we  refrain  from  stating  it  here,  in  the  hope  that 
some  alternate,  simpler,  upper  bound  improvement  might  be  obtained  later, 
in  a different  way. 
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